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Todd Springer/ Charles Gale,^'|3 Sangyong Jeon,^'|^] and Su Houng Lee^-^'j^ 

^Department of Physics, McGill University, Montreal, Quebec, H3A2T8, Canada 
^Institute of Physics and Applied Physics, Yonsei University, Seoul 120-749, Korea 
^ Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

(Dated: Jan. 17, 2011) 

A sum rule which relates a stress-energy tensor correlator to thermodynamic functions is examined 
within the context of a simple non-conformal gravity dual. Such a sum rule was previously derived 
using AdS/CFT for conformal A/" = 4 Supersymmetric Yang-Mills theory, but we show that it does 
not generalize to the non-conformal theory under consideration. We provide a generalized sum rule 
and numerically verify its validity. A useful byproduct of the calculation is the computation of the 
spectral density in a strongly coupled non-conformal theory. Qualitative features of the spectral 
densities and implications for lattice measurements of transport coefficients are discussed. 



I. INTRODUCTION 

Sum rules are powerful tools useful in the exploration 
of nonperturbative phenomena. The use of sum rules in 
this way dates back nearly three decades now [1]; but 
recently there has been some interest in applying these 
tools to the strongly coupled plasma created at the Rel- 
ativistic Heavy Ion Colhder (RHIC) [M5]. In [i[7], the 
authors used the low energy theorems of |8] to write down 
a sum rule which relates an integral over the spectral 
density to thermodynamic quantities. In general such 
a sum rule provides constraints on the spectral function, 
which could itself be used to extract transport coefhcients 
(via Kubo's formulas). Constraints or some knowledge of 
the functional form of the spectral density are generally 
needed in order to have any hope of extracting trans- 
port coefficients from the lattice. The authors of |6l [7] 
were able to argue for qualitative features of the bulk 
viscosity of the quark-gluon plasma (QGP) by combining 
their sum rule with both an ansatz for the spectral den- 
sity and lattice data. This approach was later criticized 
in [S], and later some corrections and clarifications were 
added in [TO] . This latter work also derived several other 
sum rules using Kramers-Kroenig relations. Subsequent 
works have derived additional sum rules, and examined 
the applications of such sum rules to lattice computations 
[TT| - [T5] . One such sum rule derived in jTUj was derived 
using AdS/CFT and was found to be applicable to A/" = 4 
supersymmetric Yang-Mills (SYM) theory. It is the aim 
of this paper to examine this sum rule within the context 
of a non-conformal gravity dual theory using the tools of 
the AdS/CFT correspondence [MHTT]. 

The sum rule in question is [10] 
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Here e is the energy density, and p is the spectral density 
^shcar(^) = -Im Gt'^'-'^iw). (2) 

The retarded Green's function in the "tensor" channel is 
defined as: 

G''^'''''{w) = J d^a^e™* ([T^^(a;), T^^(O)]) 6'(t) (3) 

and the subscript T = means the quantity of interest 
is evaluated in the limit of zero temperature. 

Both sides of the sum rule can be computed using 
AdS/CFT techniques. The left side depends only on 
thermodynamic quantities, which are easily evaluatable 
for the theory of interest. In order to evaluate the right 
hand side, one needs to compute the spectral density p as 
a function of w. In AdS/CFT, the differential equations 
necessary to compute spectral densities are often difficult 
to solve analytically (though in some cases analytical re- 
sults have been given in the literature jl8j). In this work 
we will solve the differential equations numerically, and 
hence our verification of the sum rule will be numerical 
in nature. 

In [10] , the authors checked that the left and right sides 
of the sum rule ([l]) are in agreement within the context 
of the (conformal) Af — A SYM theory. However, the 
authors then state that the sum rule should hold for any 
Einstein gravity dual. As shown below, this is actually 
not the case. We have evaluated the left and right sides 
of the sum rule ([T]) in a particular non-conformal gravity 
dual theory and find that the left side is not, in general, 
equal to the right side. In fact, one should expect that 
the sum rule should be corrected as 

ls-^F{s,P,v.) r ^Ap^^^^^iw). (4) 
5 n Jq w 

We will often employ the shorthand notation A to de- 
note a quantity which has the zero temperature part sub- 
tracted out. For example, 

Ap{w) = p{w) - pT=o{w). (5) 

To be consistent with currently known results, F{e, P, Vs) 
must vanish when e = 3P and = 1/3. 
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Using the same techniques as in [TU], we have been 
able to derive the correction of the left hand side of the 
sum rule in our particular non-conformal gravity dual. 
We explicitly show that the left and right sides of our 
corrected sum rule Q agree within the numerical error. 

The non-conformal theory in which we work is a sim- 
ple 5D single scalar model with an exponential poten- 
tial. This model is sometimes called the Chamblin-Reall 
model [H], and has been extensively studied in the lit- 
erature [20U26j . We emphasize that this model is not 
particularly well suited for QGP phenomenology; it has 
no conserved charge, and also has the peculiar feature of 
being both non-conformal and having a speed of sound 
which is independent of temperature. Still, we choose 
to work in this model because it is perhaps the simplest 
example of a non-conformal gravity dual where many of 
the hydrodynamic equations can be solved exactly. It is 
worth mentioning that if there is a precise field theory 
dual to this model, it is not known at present. However, 
recently it was found that the dynamics of a more compli- 
cated string theory setup (including fundamental flavors) 
were captured by an effective single scalar Chamblin- 
Reall background [27 . This may indicate a connection 
between the Chamblin-Reall background and more rigor- 
ous non-conformal deformations oi Af = A SYM theory. 

Our paper is organized as follows. In Sec. |lTj we in- 
troduce our non-conformal gravitational dual, the sin- 
gle scalar Chamblin-Reall background. In Sec. |III[ we 
present the details of the evaluation of the right hand 
side of the sum rule. This section involves introducing a 
tensor perturbation into the geometry, numerically solv- 
ing for the spectral density, and integrating the result. In 
Sec. |IV| we evaluate the left side of the sum rule using 
the known thermodynamics of the gravity background; 
it is evident that the left side does not agree with the 
right side except in the limiting case of a conformal the- 
ory. We then proceed to derive the correct form of the 
left side and present an improved sum rule where the left 
and right sides agree numerically. We conclude the paper 
in Sec. |V] In Appendix [A] we discuss the relevant sum 
rule within the context of (weakly coupled) Yang-Mills 
theory. Other technical details of our calculations and 
useful reference formulae are found in Appendices [B] - [P] 



II. GRAVITY BACKGROUND 

The theory under consideration is a 5D gravitational 
dual generated by a single scalar fieltQ 



S = 




R - - F(0) 



(6) 



^ Throughout this work, we use the "mostly plus" metric signa- 
ture. 



where k is related to the five dimensional Newton's con- 
stant, K, = SttGs. The second term is a boundary con- 
tribution, the well known Gibbons-Hawking term which 
is necessary for a well defined variational principle. The 
induced metric on the boundary is denoted by ^p,u, 
denotes the covariant derivative, and is the trace of the 
second fundamental form 



(7) 



with a unit vector normal to the boundary. 

We will assume the metric is of the "black brane" type 

ds^ = gtt{z)dt'^ + 9xx{z)dx^ + gzz{z)dz'^, (8) 

and that the coordinates can be chosen such that there 
is a black brane horizon aX z — z^ . We will often employ 
the symbol 



f{z)^^gu{z)g-^z). 



(9) 



As mentioned in the introduction, in presenting our main 
results we will specify to a particular type of metric, the 
Chamblin-Reall background. However, whenever possi- 
ble, we will keep the metric components general in hopes 
that doing so may be useful for those wishing to do anal- 
ogous calculations in different backgrounds. 

The Chamblin-Reall background can be found by as- 
suming an exponential potential of the form 



y(0) 



6 {2-5) 



■ exp 




(10) 



The potential contains a parameter 5 related to confor- 
mal symmetry breaking; the precise form of the poten- 
tial above is chosen for future convenience. The result- 
ing metric and scalar field profile which solve Einstein's 
equations are 



~f{z)de 



dz^ 

W) 



b{z) 



L 



(11) 



(12) 



1 



2(2-^) 
1-25 



(13) 
(14) 



^log[5(z)]. 

Here L is a constant which is related to the radius of 
curvature of the space, and Zh is the position of the hori- 
zon. The coordinate z runs from to Zh, with the UV 
boundary at z = 0. 

Thermodynamics and transport coeflicients have been 
studied in this setup in [501 121 123 [Ml HH] ■ The relevant 
results for our purposes are 
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1 

_ 1 
~ 3 

= 2 



2(5 



2(2 -(5) 



Ts 



(1 - 2<5) 
1 

1 

3 



(15) 
(16) 
(17) 



3 



Here we have introduced s as the entropy density, P as 
the pressure, and Vs as the speed of sound. Note that 
the parameter 5 is a measure of the conformal symmetry 
breaking; for 5 = 0, we recover the usual ^^5*5 metric, 
which is dual to a conformal field theory. One should also 
note that this setup is rather peculiar in that the speed of 
sound is constant with respect to temperature, though it 
is not necessarily equal to We will always work in 

the regime where < (5 < 1/2; in this regime the speed 
of sound is positive and less than 

III. RIGHT SIDE OF SUM RULE 

A. Tensor mode perturbations 

In order to access the two point correlation functions, 
we must add perturbations to this geometry. We assume 
a perturbation which depends only on time and the extra- 
dimensional coordinate z. 

9,1^ gfiiy + hf_,^{t,z) (18) 
<j)^^Q+d(bit,z). (19) 

In general, we might also assume a spatial dependence 
for the perturbations. Upon Fourier transform, we would 
acquire a momentum dependence of e^^'^. The sum rule 
in question involves the two point correlation functions 
at vanishing spatial momentum so we have set k to zero; 
this is equivalent to assuming the perturbation does not 
depend on the spatial coordinates x'. 

Perturbations in the 4D fluid are generally catego- 
rized into scalar, vector, and tensor modes denoting their 
transformation properties under spatial rotations. We 
will examine the tensor mode; this is the mode which 
gives access to the shear viscosity. In this case, the only 
nonzero metric perturbation is h^y, and we need not con- 
sider the fluctuation 6(j), as it does not couple to the met- 
ric perturbation in this channel. 

In order to compute the correlation functions one must 
solve the linearized Einstein equations for the perturba- 
tion's profile. Once this is accomplished, one must plug 
the result back into the action and use the prescription 
of Son and Starinets [T7] to get the correlation functions. 

The linearized Einstein equation of motion for the 
Fourier transform 

H{t, z) ^ hlit, z)=J ^H{w, z)e-™* (20) 

is, 

-J^^d.- [V^g^'H'] - w'g,.g''H = 0. (21) 



Throughout this work, we use the prime to denote deriva- 
tive with respect to the coordinate which labels the extra 
dimension (in the case at hand, z). This equation needs 
to be solved with the "incoming wave" boundary condi- 
tion which can be applied by making the ansatz 

H{z) = f{z)-^^/^Y{z), (22) 

and requiring that y is a regular function of z at the 
horizon. We have defined the customary dimensionless 
frequency 

ro = w/(27rT). (23) 

The solution for H will contain one integration constant, 
which can be related to the boundary value of H{z — >■ 0). 

Correlation functions of the operator dual to the fluc- 
tuation H can be found from the on-shell action. In 
order to access two point functions, one needs to expand 
the gravitational action to second order in perturbation 
H. Upon application of the equations of motion ("on- 
shell"), the action reduces to boundary terms, some of 
which arise due to integration by parts. We will not 
go through the details of this procedure, but they can be 
found in complete generality in P51 - I5T] . Here, we present 
the on-shell boundary terms for the theory in question: 
a gravity dual with a black brane metric supported by a 
single scalar field. 

In writing the following expressions, we have chosen to 
remove all instances of the potential with the background 
equations of motion. These equations can be found in 
Appendix |B] We have also made use of the fact that all 
backgrounds generated by scalar fields only must satisfy 

m 

gM^c^[g,,{zt^-^Pf. (24) 
/(2) 

This constraint is a consequence of the background Ein- 
stein equations. Here, ci is a constant related to the 
temperature 

" WTYgU^' ^^^^ 

Combining these expressions with the Beckenstein- 
Hawking entropy law s = gxxi^h)^^^ /"^G allows one to 
remove Newton's gravitational constant G in favor of the 
thermodynamic quantities T, s, and the constant ci. In 
computing the on-shell action, we find the constant ci 
drops out completely. The on-shell action is divergent 
and so we introduce e as a UV cutofS There are two 
contributions to the on-shell action which can be written 
in terms of the quantity 



^ One should take care to distinguish the UV cutoff e from the 



energy density e. 
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S 



Aiz),B{z) 



TsV 



dw f{z 



2n f'{z) 



y^l^A{z)H'{w, z)H{-w, z) + B{z)H{w, z)H{-w, 



(26) 



We have introduced V to indicate the spatial volume, the 
result of the integration over (Px. One contribution to 
the on-shell action is the results from the fact that bulk 
action (the first term in ([6])) reduces to a total derivative 
upon application of the equations of motion: 



^bulk — S 



-3, -Vl [gxx{z)] 



(27) 



We often employ the notation I?l to denote the logarith- 
mic derivative 



Vl[X]=X'/X. 



(28) 



In addition, there is a contribution from the Gibbons- 
Hawking term (the second term in (|6|), which is 



^GH 

In total, then 

^Total = ^bulk ^ 



S 



A,VL[mgUz)] 



(29) 



^GH — S 



l,VL[f{z)gUz)] . (30) 



In general, the on-shell action needs to be regularized 
with the addition of counter terms. However, if one is 
only interested in the imaginary part of the correlators 
(the spectral density p), this is not necessary. It is well 
known by now that the imaginary part of the on-shell 
action independent of the radial coordinate and is not 
divergent. The imaginary part i^ 



Im S- 



Total 



2i 



[5Total(w) - 5Total(-w)] (31) 



which can be written 
TsV 



Im S' 



Total 



dw f 

27T f 



H'{w,z)H{-w,z) 



- H'{-w,z)H{w,z) 



(32) 



As mentioned above, one can evaluate this quantity at 
any value of z; a convenient one is z = z^. We assume 
that our function / vanishes linearly at the horizon; 

fiz = Zh) = /o(z - zh) + 0{{z - Zhf), (33) 

with /o a constant. Furthermore, note that the incoming 
wave boundary conditions require 

H'{w, Zh)(z - Zh) = -^H{w, Zh) [1 + 0{z ~ Zh)] (34) 



^ This is because the on-shell action S is related to the retarded 
Green's function Gji, which satisfies Gji(w)* = Gji{—w). 



I 

very near the horizon. Thus, 

VTs f dw 



Im iSxctal 

H{-w,e) 



-W X 



4 7 27r 

H{w,Zh)H{-w,Zh) 
H{-w,e)H{w,e) 



H{w,e). (35) 



Here, H{w, e) is the boundary value of the perturbation 
as z — )■ 0. The prescription of Son and Starinets states 
that the two point correlation function of the operator 
dual to H is given by [11] 



- Im Gg"=^'' 
or, equivalently, 

p{w) = 



Tsxo 



H{w,Zh)H{-w,Zh) 



sw 
47r 



H{-w,e)H{w,e) 

Y{w,Zh)Y{-w,Zh) 
Y{~w,e)Y{w,e) 



(36) 



(37) 



In the case of w = 0, the only solution to the equations 
of motion which obeys the boundary conditions is H — 
constant, and thus, 



1] = lim 



p(w) 



lii-i-O w 



s 

An 



(38) 



The first equality is the Kubo formula for the shear vis- 
cosity. This is now a familiar result. 



B. Numerical computation of spectral density 

We are interested in the quantity p at finite values 
of w, and thus the equations of motion must be solved 
numerically. The first step towards this end is to pass 
to a more convenient coordinate system. We define a 
dimensionless coordinate u such that the horizon is at 
u — I and the boundary is at u = 0. One can write the 
metric as: 



a 



[-dt^fiu) + dx^ 



2(2 + 5) 

u ^-2 du^ 
a(27rr)2 J(u) 



f{u) = l-u^ 



(39) 
(40) 



The constant a is not important for our purposes (it 
drops out of the equations of motion), but for complete- 
ness, it is 



2(2 -(5) 
AttTL{1 - 25) 



(41) 



In this coordinate system, the equation of motion (21) 
becomes 



H"{u) 



1 



u{\ — v?) 



2(1 + S) 

u ^-2 



(1 



,H{u) = 0. (42) 
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Upon insertion of the incoming wave ansatz (22), the 
equation for Y is 



Y^, l + ^"(l-2^l^l, ^^, 
u (1 — u^) 



1 - ) u2fo2 

2 -r = o. (43) 



We use Mathematica's NDSolve function [35] to solve 
this equation numerically for a given value of tti and S. 
Boundary conditions must be specified. The function Y 
must be regular at the horizon in order to comply with 
the incoming wave boundary condition; the easiest way 
to apply this condition is to begin integration at some 
value of u close to the horizon and specify 



y(i-6) = i. 



(44) 



One also needs to specify the derivative Y'{u) here. Ex- 
panding the equation (43 1 in powers of (1 — u), one finds 
that the leading order term leads to the condition 



37tD^ 

^'"--'- 2(. + n,)(2-.) ^-"-» 



To summarize, our numerical method is as follows 



1. Specify a value of 5 and tr. 



2. For these values of S,Xv, use Mathematica's ND- 



Solve to numerically integrate ( 43 1 with the bound 



ary conditions (44), (45). Start the integration near 
e, and integrate down very 
e. 



the horizon at w = 1 
near the boundary at u 



3. Using the now known values of Y(e) and (37 1 one 
can determine the spectral densitjrl 



„shear 



(IT)) Ts 1 



tt) 



\Y{e)\' 



(46) 



C. Zero temperature subtraction 

The quantity that enters the sum rule is the zero tem- 
perature subtracted spectral density. One can compute 
pshcar^y^-j^^^ analytically. A gravitational metric dual to 
a zero temperature field theory is, intuitively, one without 
a black brane horizon. Returning now to our original z 
coordinates, we set f{z) — 1. In terms of b{z)'^ = gxx{z), 
the equations of motion become 

H^{z) + VLigli'izmz) + w'Hoiz) = 0. (47) 

We are using Ho{z) to denote the solution at zero tem- 
perature. We are interested in the case where g^xiz) = 
( , with n being a function of delta. 



2 

1 - 2S' 



The equation is then 



3n 
2^ 



Hl^{z) + w^Ho{z) =0. 



(48) 



(49) 



This equation is solved in terms of Bessel functions, or 
alternatively in terms of Hankel functions of the first and 
second kind: 



2 + 3n 
Z 1 



Cihi^ (wz) + C2h\^ (wz) . (50) 



(2) 



The combination (2 + 3n) /4 appears frequently in our 
calculations, and for simplicity we will use the definition 



I 



2 + 3n 



1 - 25 



(51) 



The Hankel functions of the first kind behave a,t z —>■ 
oo as '-^ e™^/-y/z, whereas the Hankel functions of the 
second kind behave as ~ e~™^' /^/z. One can think of 
the zero temperature metric as possessing a "horizon" at 
z = oo, hence we should choose C2 = so that waves are 
only traveling towards the "horizon" |17| . 

To get the correlation functions, we again need to ex- 
pand the zero temperature on-shell action to quadratic 
order in the perturbation Hq- The steps are analogous 
to those above, only now we are working in a coordi- 
nate system where —gtt — gxx — gzz- The results can be 
written in terms of the quantity 



?T=0 



A{z),B{z) 



sV f dw 
8^ .L-. 2^ 



gxxjz) 
gxx{zh) 



3/2 



[A{z)H'q{w, z)Ho{-w, z) + B{z)Hq{w, z)Ha{-w, z)] 



(52) 



* One may be worried about applying our results to this new coor- 



dinate system; in fact, the only coordinate dependent assumption 
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Finally, using the definitions for n and Z (48) and (51 1, 



The results for the bulk and Gibbons-Hawking terms are we have our final result for the zero temperature case: 



<-T=o _ 

'-'bulk — 



T=0 



<^T=0 _ <^ 

CT=0 
'-'Total 



T=0 



-s 



T=0 



3,T)l [9xx{z] 
l,VL[gUz)] 



(53) 
(54) 
(55) 



And again applying the prescription of Son and Starinets, 
we find the spectral density, 



Jxx (^* ) 



\gxx{zhj 

Hq{w, z^)Hq{-w, z^) - Hq{-w, z^)Ho{w, z*) 



3/2 



Hoiw,€)Hoi-w,e) 



, (56) 



The symbol is used to denote any particular value of 
z which we choose (again, this result is independent of 
z). In computing the finite temperature spectral density, 
we found it most convenient to evaluate the result at the 
horizon; here it is more convenient to choose z* = e. 
Let us also now specify to the case at hand with g^x = 
(L/z)«. 



Zh \ 



H'oiw, e)Hoi~w, e) - H'^{-w, e)Hoiw, e) 



Hoiw,e)Ho{-w,e) 



With the use of the solution ( 50 ) , one finds 



H'^{w,e)Ho{-w,e) - H'^{-w,e)Ho{w,e) 

_ 4|Cip(-l)-3"/^(£)'^"/^ 

Finally, one needs to employ the expansion 



(57) 



(58) 



.(1) 



r(0-t-O(a;2) + O(x'). (59) 



Recall that Z > 2 in the physical region. Putting it all 
together we find 



sw{wzhY"'^'^ 

22'+ir(02 ■ 



(60) 



Finally, one should remove Zh in favor of T. When 
the metric is written in the coordinate system (11), the 
Hawking temperature is 



T 



fjzu) 
An 



I 



2nzh 



(61) 



we have made is that the gtt component vanish Unearly near the 
horizon, and the g^z or Quu component diverge as 1/x near the 
horizon. These facts are true in both coordinate systems. 



ttTs 



2T 



(62) 



After numerically computing the spectral density at finite 
temperature, one subtracts off this piece to remove the 
UV divergences at large ir. 

One drawback of our method is that it requires a large 
degree of numerical precision; both the zero tempera- 
ture and finite temperature spectral functions diverge at 
large tn, but their difference approaches zero. One needs 
to compute both spectral functions to a high degree of 
precision before performing the subtraction to get the de- 
sired result. In the future it would be desirable to build 
this subtraction into the numerics, to avoid this need for 
very high precision numerics. 



D. Numerical results for spectral density 

We have numerically computed the (zero-temperature 
subtracted) shear spectral function for 6 = 0, 0.1, 0.2 and 
0.3. Some sample results are shown in Fig. [ij We have 
computed the spectral functions out to a large value of 
Id; we cease our numerical computation when the oscil- 
lations have reached 0.1% of their maximum value. The 
qualitative behavior of this spectral function is quite in- 
teresting. The spectral density oscillates around the zero 
temperature result, with the oscillations eventually dying 
out as one moves to higher frequencies. This oscillation 
phenomenon has been noticed before in previous com- 
putations of the spectral density in A/" = 4 SYM theory 
[551 ED ■ These damped oscillations are thought to be a 
reflection of the pole structure of the retarded correlation 
functions in the complex plane |18l 135) . 

When moving to a non-conformal theory, we notice a 
qualitatively new behavior. As one increases the non- 
conformal parameter 5, the oscillations become initially 
larger and larger, and the spectral function takes a longer 
time to settle down to its zero temperature value. This 
is especially evident for the case of 5 = 0.3 where the 
oscillations are so large that it does not fit nicely on a 
plot with those shown in Fig. [T] We are unsure of the 
physical interpretation of this behavior. However, we 
can use our spectral functions to compute Euclidean cor- 
relation functions, quantities that are computed on the 
lattice (see for example [SHI 133 )■ The relation between 
the spectral function and the Euclidean correlators is 



Ge{t) 



dwp{w) 



cosh [w{t - /3/2)] 
sinh [w/3/2] 



(63) 



where r is the Euclidean time variable, which has period 
/3 = 1/T. Using our numerical results for the spectral 
density, we plot the results for the Euclidean correlation 
functions in Fig. [3j It is interesting to note that the 
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Zero temperature subtracted shear spectral density 




Euclidean Correlation Functions 




0.3 



0.4 



0.5 
tT 



0.6 



0.7 



FIG. 1. Plots of the zero-temperature subtracted spectral 
density versus frequency for several values of 5. For large w, 
the spectral density always approaches the zero temperature 
result. For small w, the slope of the spectral density always 
approaches the same result which is confirmation that the 
shear viscosity takes on the universal value r]/s = l/4n. 

Zero temperature subtracted shear spectral density 




2 3 4 

w/(2 It T) 

FIG. 2. Plot of the zero-temperature subtracted spectral 
density versus frequency for 5 = 0.3. One should note the 
qualitative differences in the plots as one increases S. For 
larger values of S the oscillations grow larger and more fre- 
quent, and take a longer time to die off. Note the difference 
in the axis scaling between this plot and Fig. [l] 



value of S (which is proportional to the bulk viscosity C) 
affects the magnitude of the Euclidean correlation func- 
tion quite strongly. This may have implications for lattice 
measurements, since the Euclidean correlation function is 
directly measured there. The gravity dual in which we 
are working is only toy model of a non-conformal theory, 
and hence we will not make any attempt to extract a 
value of the bulk viscosity from the lattice data. How- 
ever, if the qualitative behavior noted here persists in 
more realistic holographic models of physical gauge the- 
ories, it suggests that one could possibly gain insight into 
the value of the bulk viscosity from the tensor correlation 
function considered here. This could prove to be quite 
useful, since it would provide an independent measure- 
ment of the bulk viscosity which is usually extracted via 
the Kubo relations in the bulk channel. 

To evaluate the sum rule, one must integrate the spec- 
tral functions. Technically, we use Mathematica to per- 



FIG. 3. Plot of the Euclidean correlation function associated 
with the shear spectral density as a function of the Euclidean 
time T for various values of S. The value of 5 (and hence the 
value of the bulk viscosity in our model) has a strong effect 
on the shape and magnitude of these functions. 



s 


(LHS of Sum rule) x ^ 


(RHS of Sum rule)x^ 


0.0 


^ ^ 0.300... 


0.300... 


0.1 


^ ^ 0.318... 


0.326... 


0.2 


i « 0.333... 


0.357... 


0.3 


^ ^ 0.353... 


0.395... 



TABLE I. Results for the left and right sides of the sum rule 
given in [TD]. The left side values come from ([T]) and (15 1. 
The right side values come from numerically integrating our 
results for the spectral density. 



form a cubic interpolation between the points which are 
spaced at intervals of tr = 0.01, and numerically integrate 
the resulting function. The results are given in Table |l] 



IV. LEFT SIDE OF SUM RULE 

If we take the sum rule ([I]) at face value, we can easily 
evaluate the left hand side using (15). The results are 
shown in Table |l] Clearly, there is substantial disagree- 
ment between the left and right sides. Notice that in the 
case of (5 = 0, the two sides are in agreement, which was 
also the conclusion of [TD] . However once we deviate from 
conformality, differences appear. The error for S — 0.3 is 
greater than 10% which is not accounted for by numerical 
error. 

Thus, we come to the conclusion that the left hand 
side of the sum rule must be modified. Following jlOj, 
the sum rule is, more generally 



{w — ioo) 



[w = 0) 

OO 7 

w 



{w). (64) 



The left side of (64) can be computed using AdS/CFT 
techniques. The general method we use is the same as 



8 



[lOj . except we are now working in a more general back- 
ground, and hence the equations of motion and action 
are modified. 

We win first compute the term AG'|J'°^'' = ioo). Let 
us define 



(65) 



fn terms of the dimensionless variable q = «tt), the limit 
is q — > oo. Note that this limit can also be achieved by 
taking T — > 0, or equivalently — )■ oo. This, in turn is 
equivalent to examining the metric in the regime near the 
boundary z — 0. Our first step is to transform the met- 
ric to Fefferman-Graham like coordinates; the relevant 
details of this transformation can be found in Appendix 
[Cl Near the boundary, the metric can be written as 



-1 



3-4(5 



4(5 



4-4(5 

21 



(66) 



One must now solve the equation of motion using this 
metric. Since we are working in the regime of Q — ^ oo, we 
will denote the relevant fluctuation as . Solving the 
equations of motion in this regime is rather technical, but 
details can be found in Appendix [Dj The near boundary 
solution is 



(67) 



KiiQz) - - 

Zh 



2 

Qzh 



3r(0 

4(5 - 4^) 



O 



Here C3 is a normalization constant, and Ki is a Bessel 
function. The reader is referred to Appendix [D| for more 
details. 

To compute the correlator, we use the results for 
the on-shell action given previously ( 30 ) , the Son and 
Starinets prescription gives 



Gr{w = ioo) ~ Ts lim 



f{e) H°-'iQ,e) 



of'ie) i/-(g,e) 



+ G'^'^(w = ioo). 

(68) 

Here, the second term with the superscript CT denotes 
the contact terms which arise from the part of the on-shell 
action which contains no derivatives of H. We will have 
more to say about this term in a moment. Our solution 
(67) for H°° contains two terms; the first term is the 
result if one takes the horizon — 00; in other words, 
it is the zero temperature result. Upon subtracting the 
zero temperature piece, we have 



AGr{w 



— Ts lim 



/(e) H^'{Q,e) 



o/'(e) H^{Q,e) 
- ioo). 



(69) 



Noting that 



Zh 



with ^ a constant, and 



z2l / ~4l 



(70) 



(71) 



we have. 



AGfi{w ~ ioo) = Ts- 



AG'i'^{w = ioo). (72) 



Finally, with the fact that 



<(Q,e)^l(^^) m+Oie), (73) 



we come to the result 



AGr{'w = ioo) = Ts 



2(5 - 4(5) 



AGg^(u' = ioo). 



(74) 

Ostensibly, we are not finished as the left side of the sum 
rule also contains a term AGb{w = 0), and we also have 
to deal with the contact terms arising from the part of 
the action containing no derivatives of H. There is also 
the issue of counter terms which could be added to reg- 
ularize the action. The claim of [10] is that the contact 
terms above precisely cancel the contribution at w = 0, 
and that all counter terms will cancel if the zero tem- 
perature subtraction is done properly. We will not prove 
this claim here, but the agreement between our derived 
formula and our numerical results is an empirical justifi- 
cation for this claim. We believe that showing this explic- 
itly should be mostly straightforward given our solutions 
and our expansions for the on-shell action; however the 
zero temperature subtraction can sometimes be a sub- 
tle issue (see for example, [3S]). The final result of our 
analysis is 



AGr{w = ioo) - AGn{w = 0) = Ts 



2(5 - 4:6) 



(75) 



Using now our improved formula we compare the left and 
right sides of the sum rule in Table |lTj As one can see, 
the agreement is much better, and the numerical results 
agree with the analytic ones to at least three significant 
figures. The fact that the error increases with 6 is not 
surprising, since the numerics become more challenging 
as (5 increases. This is because the power divergences 
of p and pt=o are worse, and hence greater numerical 
precision is required. Furthermore, as S increases, the 
oscillations take longer to die out, and hence one must 
compute Ap to a larger value of ro in order to achieve the 
same accuracy. 

The largest source of error in these computations is due 
to the fact that we only compute p up to a value rDmaxj 
and our numerical integration stops at this value. As 
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5 


{New Sum rale LHS)x 


(Sum rale RHS) x 


0.0 


^ ^ 0.300000... 


0.3000(04)... 


0.1 


if « 0.326087... 


0.3260(86)... 


0.2 


^ « 0.357143... 


0.3571(28)... 


0.3 


If « 0.394737... 


0.395(036)... 



TABLE II. Comparison of the left and right sides of the im- 
proved sum rule (771 which is now corrected to include the 



non-conformal effects of the theory. The left side is computed 
from ( 75 1 , while the right side is computed by the numerical 



integration over the spectral density. The parentheses denote 
an estimate of numerical error which is explained more com- 
pletely in the text. 



shown in [TU], the zero temperature subtracted spectral 
density is an oscillating function which dies out nearly 
exponentially at large fo. We are neglecting the integral 
of this function from tn — fOmax to It) = oo. One can 
estimate the contribution of this "tail" by fitting the last 
few oscillations near ttimax to a damped sine curv £0 of the 
form 



tX) 



-6ro 



sin (clD -t- d) . 



(76) 



Once the parameters a,b,c and d are found from the fit, 
one can then integrate this function from fOmax to oo 
to estimate its contribution. We find that the order of 
magnitude of this tail contribution is ^ 10^^ for 5 — 
0,0.1, and 0.2. The contribution of the tail appears to 
increase with 5, and is roughly ~ 10"'' for 5 = 0.3. This 
is the significance of the parenthesis in the table of our 
numerical results; it is expected that the contribution of 
the tail will affect the numbers inside the parenthesis. 
However, it is quite clear that our analytic results agree 
with the numerical ones within our estimated numerical 
error. 

The sum rule for the particular non-conformal theory 
in which we are working is 



Ts 



dw 



(77) 



2(5- 4(5) _ 

It is desirable to write the left side in terms of thermo- 
dynamic observables, in hopes that our sum rule could 
be applicable to other theories beyond those considered 
here. Unfortunately, there is not a unique way of doing 
that in our case, as one could write either 



S = 



3P 



2e 



or S=l(l-3v^^ 



(78) 



^ This functional form was chosen for the purpose of error estima- 
tion only, and appears to fit the numerics quite well. In fact, 
the form we have chosen matches an analytical calculation of a 
current-current correlation function in 1181 . However, we stress 
that the exact functional form of Ap'*'^"'"' is not known analyti- 
cally at large vo, though it should be possible to calculate it using 
methods similar to those given in this paper. 



This is a consequence of the fact that the speed of sound 
is independent of temperature in this model. One could 
perhaps gain more information by investigating this sum 
rule in a gravitational dual theory where the speed of 
sound depends on temperature (e.g. 1551 15^ ) 

We elect to write the left side in terms of e and P only. 
Then, the sum rule can be written 



eie ■ 



P) 



2{e + 2P) 



2 

./O 



dw 



(79) 



This sum rule reduces to ([T]) in the case of e = 3P. In 
other words, we have found that the function F{e,P,Vs) 
defined in ^ can be written 



F(e,Pv,, 



eie - 3P) 
10(£ + 2P)' 



(80) 



V. CONCLUSION 

In this work, we have examined a sum rule involving 
a particular two point function of the energy-momentum 
tensor. A version of this sum rule was derived in [lOj . 
but we have explicitly shown that the sum rule given 
there is not valid for all Einstein gravity dual theories. 
We have provided a non-conformal generalization of this 
sum rule; the main result of this work is (79). We have 



numerically verified that the left side of our improved 
sum rule equals the right side with an accuracy greater 
than 0.1%. Whether our result is applicable to other non- 
conformal gravity duals, or even to non-conformal field 
theories should certainly be tested. To this effect, we 
have examined the sum rule in Yang-Mills theory in Ap- 
pendix |A] In addition, it is also important to numerically 
verify other sum rules given in |10] . These investigations 
are currently underway. 

Finally, in the course of our investigation of the sum 
rule, we computed the spectral density in the tensor 
channel at various values of our non-conformal deforma- 
tion parameter S. The behavior of the spectral density 
and the associated Euclidean correlation functions ex- 
hibit interesting qualitative behavior as a function of 6 
as shown in Figs. [l]-|3] While the qualitative features of 
the spectral density and Euclidean correlation functions 
were found in the conformal case of A/" = 4 SYM theory 
[551 15i] , to the best of our knowledge the change in these 
functions as a result of non-conformality has not been 
examined before. In particular, our results seem to sug- 
gest that the bulk viscosity (which is proportional to S) 
may have strong effects on the shape and magnitude of 
the tensor Euclidean correlators measured on the lattice. 
This is noteworthy because the bulk viscosity is usually 
found from the Kubo relation which involves a different 
correlation function. It remains to be seen whether these 
qualitative features will persist in more detailed gravita- 
tional dual models which attempt to capture more fea- 
tures of the quark-gluon plasma. 
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Appendix A: Sum rule in Yang-Mills theory 

We will use this section to make some comments on the 
shear sum rule in (weakly coupled) Yang-Mills theory. A 
tentative form for this sum rule was written in [TU^; but 
here we will perform the calculation in a slightly different 
way. 

The left side of the sum rule contains correlation func- 
tions at large frequency. Such quantities can be calcu- 
lated using the operator product expansion (OPE). In 
the limit of w — >■ c», one only needs to know the OPE 
to leading order. In [JDIj the leading order OPE was cal- 
culated for the correlator of interest, with the result (in 
Euclidean signature) 

G^h'=-(g ^ oo) = ^ (g^ _ (t44) + I (F') (Al) 

Here, F^^ is the gluon field strength tensor, and = 
P^u^a'^ ■ The stress-energy tensor is denoted by T'"^, 



rp pa pa A _ p2 



(A2) 



and the index "4" pertains to (74 = iw. In the limit 
of zero spatial momentum , and returning to Minkowski 
signature, we find 

A D 

All that remains is to express these quantities in terms 
of thermodynamic ones; the arguments here follow [?T] . 
With the assumption of a perfect fluid, and with the use 
of the trace anomaly we can write 



(TO") 



(A4) 



(^^') = -^W = J^(^-3P), (A5) 



/3(.9) 



where the scale dependence on the left hand side of ( A5 ) 



is transformed to the scale dependence of the coupling 
constant. For notational convenience, we define 



shear 



(A6) 



This is the quantity that will enter the left side of the 
sum rule for pure gluodynamics, the low energy term 
AGbXw = 0) vanishes as shown in Appendix A of jlO| . 
Note that the second term in ( A8 1 which is propor- 
tional to aj^ was not given in [10 , though it appears 
to have been independently noticed in a recent paper 
|13) . This additional term originates from the (T^) 
In |— 



term m 



Al 



this term is argued to be a con 
tact term due to the fact that it renders the correlation 
function non-transverse. If this is the case, this addi- 
tional term should not be present in the sum rule, since 
all contact terms will cancel out due to the subtraction 
AGr{w — ioo) — AGi?(w = 0). Despite the claim made 
in [4G| , we are unaware of any explicit calculation which 
shows that AG^'"'''{w = 0) ^ e - 3P. For this reason, 
we currently choose to include the additional term in the 
sum rule, though it is clear that more work should be 
done to address this issue in the future. 

In perturbation theory, e — 3P ~ ©(a^), so in total, 
the second term on the right side of (A8) is 0{as). Let 
us make a few comments on what may happen to this 
expression as we go beyond leading order in perturbation 
theory. We expect that the Wilson coefficients in (Al) 
will, in general, contain corrections of 0{a 
words, we could expect that (Al) becomes 



In other 



G^hcar^^ ^ oo) = — (1 + ai«.) {qj - q^) (T^^) 



+ -(l + a2«.)(P') 



(A9) 



with constants ai and a2 to be determined from a one 
loop calculation. Following through the previous argu- 
ments, we would arrive at the expression 



- GS°" = ^ (1 + «i«s) s + (e ~3P) + (al) . 

(AlO) 

or, equivalently. 



G: 



ihcar 



1 



{I + aias) (e + P) 



(All) 



47r 



(e -3P)+0 {al 



We find 



where again we have used the fact that e — 3P is O (a^) . 

We find it interesting that the left side of the sum rule 
in our strongly coupled model has a different dependence 
on £ and P than the left side of the sum rule in weakly 
coupled Yang-Mills theory. It is not clear whether this 
has any implications for the possibility of a gravity dual 
of weakly coupled Yang-Mills theory. 



'shcar 



3^ 3/3(g) 



(£ - 3P) . 



To leading order, this becomes: 

2 4 
3^ ^ im^as 



^shear _^ 



(A7) 



2 An 

' ie-3P)- (A8) 



Appendix B: Background Equations 

For the black brane type metric generated by a single 
scalar field, with a z denoting the radial coordinate, the 
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background equations of motion can be written 

-0 (Bl) 

3 



2^ 



1 



'9 



dV 



cl>=cl>o(z) 



(B2) 

(B3) 
(B4) 



where / and are defined as in the text (|9| and (28) 
respectively. 



Appendix C: Transformation to FefFerman-Graham 
like coordinates 



In this section, we detail the transformation of the 
Chamblin-Reall metric to a Fefferman- Graham like coor- 
dinate system. Fefferman-Graham coordinates are useful 
in many respects for asymptotically anti de-Sitter met- 
rics; these coordinates are defined so that the metric 
takes the form 

dp) , , 
^, (CI) 



where the indices /i and v run over the four coordinates 

t,x\ 

Note that the Chamblin-Reall metric is not asymptot- 
ically anti de-Sitter (except for the case of 5 = 0), hence 
we will make a slight generalization of the Fefferman- 
Graham coordinates, where the metric is written 



ds' = 



{g^,^dx^'dx^' + dP) 



(C2) 



It is not always possible to solve for (jf^^ analytically, but 
one can find its behavior near the boundary as an expan- 
sion in the radial coordinate. To transform our metric 
to this form, we begin with the Chamblin-Reall metric 
written in the coordinate system [Tl] and apply the coor- 
dinate transformation 1421. 



dz^ 



^) dz' 

z 



(C3) 



Enforcing the fact that near the boundary, z = z, we can 
integrate this equation to find 

dz _ dz 

For the time being, let us assume that n ^ 2, (though 
we will be able to relax this assumption later on) . Then, 
performing the integral we find 

dz 



n 

1 

2 



A- 9 



7/2 



(C5) 



dz. 



(C6) 



We have taken the limit e — freely since the integral 
converges. Near the boundary, one can expand the inte- 
grand to fincj^ 



z^-t <^ 1 + ( 1 



O 



n\ 1 



2/ Al -n + 2 \Zh 



(C7) 



(Note that f{z) = 1 — {z/zh) .) Then, to this order we 
have, 



41- n + 2 \ z/i 



1 / Z \ z^'- 



41- n + 2 \Zh 



O 



.Al 



In terms of (5, this is. 



z«zO- 



1 - 25 




(C8) 



(C9) 



(CIO) 



Strictly speaking, we derived this result assuming n ^ 2, 
or equivalently (5^0. However, if one redoes the analy- 
sis for (5 = 0, one finds precisely the same result as that 
given by (CIO), thus (CIO) is valid for all values of 5 in 
the physical regime < 5 < | . Applying this coordinate 
transform, we find the Fefferman-Graham like represen- 
tatio n of the Chamblin-Reall metric which is given in the 
text (661 to order ^^(z^'/^^O- Since we are only interested 
in the near-boundary dynamics, we need not worry about 
higher order terms. 



Appendix D: Calculation of H°° 

In this section we present some of the technical de- 
tails involved in solving the equations of motion in the 
large w regime. Throughout this section we work in the 
Fefferman-Graham coordinate system, which is explained 
in detail in Appendix [C] Perhaps the easiest way to get 
the relevant equation is to introduce a scaling parameter 



21 



2/ 



H^{i) + \Hr{z) 



making the ansatz 
(Dl) 



A, by replacing (^j^ 



and expanding the equation of motion ( |2l[ ) in A. The 
lowest order equation gives the T = equation for Hq 
(with w"^ replaced by —Q^)- 



(Ho 



S{1 - 26) 



(H^y - Q^H^ = 0. (D2) 



We assumed 21 — ^ ^ — 1. This is valid provided 5 5/2 which 
is outside the physical regime of interest 
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The first order equation is 
3 



z(l - 26) 

21 

Z ' 
Zh 



1 



(D3) 



1 - 5 



The solution for Hn is 



H^iQ,u) = Csu'KiiQuzh) + C^u'IiiQuzh), (D4) 

where we have defined ii = z/zh for convenience. /; and 
Ki are the modified Bessel functions of the first and sec- 
ond kind respectively. Regularity at z — oo [101 E] 
requires C4 = 0. 



The homogeneous part of the first order equation ( D4 ) 
is the same as the zero-temperature equation, hence this 
equation can be solved with the use of a Green's function. 
Defining the two homogeneous solutions as 



yi{u) = Cau^KiiQuzh) 
2/2 (") = Czv!- Ii{Quzh) , 



(D5) 
(D6) 



The solution is 



CI 



vi{u) / V2{t)g{t)dt 



+ y2{u) / yi{t)g{t)dt\ 



(D7) 



where 



u 



1-5 

u 
1-5 



4 

3-4(5 



{Qzhfyi{u)-Sl 



u 



(D8) 



In order to compute the correlation function, we will only 
need the leading term at the boundary -u = 0, 



CI 



y2{u^Q) 1^ yi{t)g{t)dt} (D9) 



which can be written: 



Hi{u 0) 



1 



CsiS-l) 



u^Qzh 



IV{1) 

dt \. (DIO) 



The first term in the integrand requires the integral 

and the second term in the integrand is a total derivative 
which reduces to the boundary term ^ ?/i(0)^, (note that 
2/1 (u cxi) approaches zero exponentially). In total, the 
result is 



C^ f ~z 



;2 



«~(=^-'')-4(l-.)V20.J 

3-4(5r(l + 20\/^ <5M' 



4 r(i + 



(D12) 



The identity 



r(;-Ki)rN 



_ 0f 1 

2 ) ' 

can be used to simplify the result to: 



;2 



2z 

04 



r{2l + 2) (D13) 



(D14) 



4(5 - 4(5) 



In summary, the solution for H for large values of Q, 

H{w = ioo) = Hao can be written 



(D15) 



i?°°(z)«C3 1 - I X 



Ki{Qz) 



z \ W 2 y 3r(/) 



^ZhJ \QzhJ 4(5-4(5) 
This is the solution which is quoted in the text. 
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